Collective Dynamics of Lipid Membranes studied by Inelastic Neutron Scattering 
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We have studied the collective short wavelength dynamics in deuterated DMPC bilayers by in- 
elastic neutron scattering. The corresponding dispersion relation Tuj(Q) is presented for the gel and 
fluid phase of this model system. The temperature dependence of the inelastic excitations indicates 
a phase coexistence between the two phases over a broad range and leads to a different assignment 
of excitations than that reported in a preceding inelastic x-ray scattering study [Phys. Rev. Lett. 
86, 740 (200f)]. As a consequence, we find that the minimum in the dispersion relation is actually 
deeper in the gel than in the fluid phase. Finally, we can clearly identify an additional non-dispersive 
(optical) mode predicted by Molecular Dynamics (MD) simulations [Phys. Rev. Lett. 87, 238f0f 
(2001)]. 

PACS numbers: 87.14.Cc, 87.16.Dg, 83.85.Hf, 83.10.Mj 



The collective dynamics of lipid molecules are believed 
to affect significantly the physical properties of phospho- 
lipid membranes [l|, |2| ■ In the context of more complex 
biological membranes, collective molecular motions may 
play a significant role for different biological functions. 
For example, correlated molecular motions of the lipid 
acyl chains and the corresponding density fluctuations 
in the plane of the bilayer are believed to play an im- 
portant role for the transport of small molecules through 
the bilayer 0. Molecular vibrations, conformational dy- 
namics and 'one particle' diffusion in the plane of the 
bilayer can be studied by a number of different spec- 
troscopic techniques covering a range of different time 
scales such as incoherent inelastic neutron scattering Q| 
or nuclear magnetic resonance Q . Contrarily, few exper- 
imental techniques are able to elucidate the short range 
collective motions mentioned above. To this end, Chen 
et al. have recently presented a seminal inelastic x-ray 
scattering (IXS) study 0, of the dispersion relation 
hLu(Q r ) which quantifies the collective motion of the lipid 
acyl chains as a function of the lateral momentum trans- 
fer Q r . The basic scenario is the following: At small 
Q>, longitudinal sound waves in the plane of the bilayer 
are probed and give rise to a linear increase of u> ex Q r , 
saturating at some maximum value ('maxon'), before a 
pronounced minimum f2g ('roton') is observed at Qo ~ 
1.4 A" 1 , the first maximum in the static structure factor 
S(Q r ) (the inter-chain correlation peak). Qualitatively, 
this can be understood if Qo is interpreted as the quasi- 
Brillouin zone of a two-dimensional liquid. Collective 
modes with a wavelength of the average nearest neigh- 
bor distance 27r/Qo are energetically favorable, leading 
to the minimum. At Q r values well above the minimum, 
the dispersion relation is dominated by single particle 
behavior. A quantitative theory which predicts the ab- 
solute energy values of 'maxon' and 'roton' on the basis 
of molecular parameters is absent so far. However, the 



dispersion relation can be extracted from Molecular Dy- 
namics (MD) simulations by temporal and spatial Fourier 
transforming the molecular real space coordinates [8j. 

Here we present a first time inelastic neutron scatter- 
ing (INS) experiment on the collective dynamics of the 
lipid acyl chains in the model system DMPC -d54. Note 
that by selective deuteration of the chains, the respective 
motions are strongly enhanced over other contributions 
to the inelastic scattering cross section. We have mea- 
sured the dynamical structure factor S(Q r ,w) in the gel 
(L^) and fluid (L Q ) phases, and have investigated the 
temperature dependence of the excitations in the disper- 
sion minimum in the vicinity of the main phase transi- 
tion [19| . While qualitatively we obtain similar curves as 
Chen et al. [6J, several results are significantly different 
and in striking contrast to the earlier study. Notably, the 
minimum in the dispersion relation measured by INS is 
deeper in the gel than in the fluid phase flo.Lfl < fio,La, 
while the opposite is reported in [6j. Secondly, the INS 
data give much smaller peak widths of around 1 meV, 
indicating a much smaller damping, while Chen et al. re- 
port values in the range of Aoj s = 4 — 6 meV. Finally, we 
observe a second non-dispersive mode at higher energy 
transfer of 14 meV, in addition to the dispersive sound 
mode. This mode was not observed by Chen et al., pos- 
sibly due to the low signal-to-noise ratio in IXS at high 
energy transfer. It has been predicted by a recent Molec- 
ular Dynamics (MD) simulation and can be attributed 
to the terminal (methyl) carbons of the chains 8]. The 
measurements were carried out on the cold triple-axis 
spectrometer INI 2 and the thermal spectrometer IN3 at 
the high flux reactor of the Institut Laue-Langevin (ILL) 
in Grenoble, France. The main differences with respect 
to IXS are related to the energy-momentum relation of 
the neutron versus the photon probe, strongly affecting 
energy resolution, accessible (Q,w) range |2fJ and the 
signal-to-noise ratio. In the present case, the energy of 
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the incident neutrons was in the range of the excitations 
(some meV) resulting in a high energy resolution up to 
^300 /ieV, in comparison to 1.5 meV of the IXS exper- 
iment. A better energy resolution in combination with 
a smaller ratio between central peak and Brillouin am- 
plitudes leads to very pronounced satellites, see Fig. Q 
This is of particular advantage for the identification of 
peaks in the central part of the dispersion relation, as 
well as for the experimental verification of the predicted 
non-dispersive mode at high energies. 

Deuterated DMPC -d54 (deuterated 1,2-dimyristoyl- 
sn-glycero-3-phoshatidylcholine) was obtained from 
Avanti Polar Lipids. Highly oriented membrane stacks 
were prepared by spreading a solution of typically 35 
mg/ml lipid in TFE/chloroform (1:1) on 2" silicon wafers, 
followed by subsequent drying in vacuum and hydration 
from D2O vapor |9j- Nine such wafers separated by small 
air gaps were combined and aligned with respect to each 
other to create a 'sandwich sample' consisting of several 
thousands of highly oriented lipid bilayers (total mosaic- 
ity of about 0.6 deg), with a total mass of about 300 
mg of deuterated DMPC. The samples were kept in a 
closed temperature and humidity controlled Aluminum 
chamber. Hydration of the lipid membranes was achieved 
by separately adjusting two heating baths (Haake, Ger- 
many), connected to the sample chamber and to a heavy 
water reservoir, hydrating the sample from the vapor 
phase. Temperature and humidity sensors were installed 
close to the sample. Before the measurements, the sam- 
ples were equilibrated at T=55 °C for about two hours 
to anneal defects. The scattering vector, Q, was placed 
in the plane of the membranes (see Fig. (a) for a sketch 
of the scattering geometry) to measure the static struc- 
ture factor S(Q>) in the plane of the membranes as well 
as the dynamic structure factor S(Q>, ui) in the same run 
without changing the setup. 

A typical energy scan (raw data) is shown in Fig. ^ 
(b), taken at T=20 °C, in the gel phase of the bilayer 
at Q r =1.0 A -1 (all data are corrected for temperature 
by the Bose factor). The inset shows the excitations of 
the bilayer in the gel and the fluid phase (at T=30 °C, 
nine degrees above the phase transition temperature [il| 
). Position and width can easily be determined from 
these well pronounced peaks. The inelastic scans can be 
evaluated by the generalized three effective eigenmode 
theory (GTEE) @ 0, El, 



using the following function 



for least-square fitting 

S(Q,u) = 1 ( A _ 
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The model consists of a heat mode, centered at u> = 
meV (Lorentzian with a width T^), and two sound 
modes, represented by Lorentzians at uj = ±uj s and 
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FIG. 1: (a) Schematic of the scattering geometry, (b) En- 
ergy scan in the gel phase of DMPC (T=20 °C), measured 
at Q r =1.0 A -1 (IN12 data). The scattering is composed of 
the central elastic peak, the broad quasi-elastic background, 
and symmetric satellites indicating the excitation. The in- 
set shows a zoom of the satellite peaks in both gel and fluid 
phase (T=30 °C). (c) Scan over a higher energy transfers 
range (Q>=3.0 A -1 , IN3 data), showing both the dispersive 
excitation and the non-dispersive (optical) excitation. 



damping 7 S [t!|10|. From width of the central mode, and 
width and position of the Brillouin lines, the thermal dif- 
fusivity, the sound frequency, and the sound damping and 
can be determined, respectively, within the framework of 
a hydrodynamic theory. To fit the neutron data, we have 
to add an additional Lorentzian component describing 
the broad quasielastic contribution associated with in- 
tramolecular degrees of freedom and incoherent scatter- 
ing, not seen by IXS. The solid line in Fig. ^ (b) is a 
fit to the GTEE model with an additional Lorentzian 
component. Figure ^ (c) shows an energy scan in the 
gel phase at Q r =3.0 A" 1 , up to an energy transfer of 30 
meV. Aside from the dispersive excitation due to in-plane 
density waves, a second non-dispersive (optical) mode is 
observed at about lj=14 meV with a width (FWHM) of 
about 13 meV, corresponding exactly to the predictions 
by Tarek et al. |8(, and attributed to the methyl ends of 
the acyl chains. 

The dispersion relation in the gel and the fluid phase is 
shown in Fig. [21(a), after evaluating the energy positions 
of the satellite (Brillouin) peaks measured in a scries of 
constant Q-scans in the range from Q=0.7 A -1 to Q=3.0 
A~h The fluid dispersion has been measured at T=30 
°C, far above the phase transition and the regime of so- 
called anomalous swelling in the fluid phase of the DMPC 
bilayers. Fig. [3 (b) shows energy scans in the dispersion 
minimum, measured with an enhanced energy resolution. 
These curves cannot be fitted by only one excitation. 
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FIG. 2: (a) Dispersion relations in the gel and the fluid phase 
of the DMPC bilayer. (b) Energy scans in the dispersion 
minimum at Q=1.5 A -1 for temperatures T=20, 23 and 30 
°C . 

Rather two contributions have to be assumed, which we 
identify from their respective T-dependence, as discussed 
below. As a consequence, the dispersion relation in the 
minimum goes down to smaller values in the gel than in 
the fluid phase, distinctly different from [6j. While the 
'maxon' and the high-Q range are energetically higher in 
the gel than in the fluid phase (due to stiffer coupling 
between the lipid chains in all-trans configuration), Oo, 
the energy value in the dispersion minimum, is actually 
smaller in the gel phase, roughly analogous to soft modes 
in crystals. 

To investigate the temperature dependence of the col- 
lective dynamics, we carried out simultaneously inelastic 
and elastic measurements in the most interesting range 
around the minimum in the dispersion relation and cor- 
respondingly the maximum of the (elastic) chain corre- 
lation. Temperature was varied between T=20 and 40 
°C, to cover the range of anomalous pseudo-critical be- 
havior. Although the phase transition is of first order 
[TlIlT^ . significant (continuous) changes (anomalies) are 
observed close to the phase transition in many thermo- 
dynamic quantities, e.g. in the equilibrium distance d 
(so-called anomalous swelling). The reason for this phe- 
nomenon is still a matter of debate. To illustrate the 
characteristic changes in the (static) chain correlation 
peak, elastic scans at temperatures of T=20 , 23 and 
30 °C are plotted in Fig. |3|(a). Upon heating from the 
gel to the fluid phase, the peak position shifts to smaller 
Q r -values (larger average next neighbor distances) and 
the peak broadens, indicating a less well ordered packing 
of the acyl chains |l3j , which is quantified in Fig. [3] (b) in 
terms of fitted peak position Qo(T), peak intensity and 
correlation length £ r (T)=l/HWHM. 

Over the same T-range, Fig. |3] (b) shows representa- 
tive energy scans taken at constant Q r = 1.5 A -1 , in 
the dispersion minimum. In the gel phase at T=20 °C 
we find an excitation at u> s =-l meV, which we associate 
with the gel phase. At higher temperatures, the gel exci- 




FIG. 3: (a): Elastic structure factor showing the inter-chain 
correlation peak, measured at temperatures T=20, 23 and 30 
°C. (b): Peak position, intensity and correlation length as 
fitted for all measured temperatures. 



tations decreases in intensity, while a second excitation at 
energy values of lo s =-2 meV is found to increase with T 
and can thus be attributed to the fluid phase. As the ex- 
citations are symmetric around the central peak we are 
sure not to deal with spurious effects. The assignment 
of the excitations to the particular phases is justified by 
their temperature dependence: In each phase there is 
a dominant excitation. Both modes are clearly disper- 
sive and change energy position when moving out of the 
minimum (see Fig. HJ(a)). We find small traces of the 
'fluid excitation' already in the gel phase. At T=30 °C, 
far in the fluid phase, there is still the 'gel excitation', 
indicating a coexistence of fluid and gel domains. This 
coexistence is only observed in the range of the disper- 
sion minimum which coincides with the maximum of the 
static structure factor. Note that at the same time, the 
elastic scans do not show coexistence of two phases, but 
the typical well known behavior: While the transition is 
of first order, a pseudo-critical swelling, i.e. a continuous 
change of the interlamellar distance above T c , is observed 
for DMPC and other lipids [U H El] The changes 
in d are accompanied by corresponding changes in the 
mean distance between chains, see Fig. 01(b). For both 
quantities, the elastic diffraction data show a continuous 
change, while the present interpretation of the inelastic 
curves indicates a coexistence of gel and fluid domains, in 
particular also far above the transition in the fluid phase. 
A crucial point is the size of these domains: The coexis- 
tence of macroscopic domains with sizes larger than the 
coherence length of the neutrons in the sample would 
lead to a peak splitting of the chain peak. Therefore, 
the domain sizes must be smaller than a few hundred A 
estimated for £„. Fluid domains in the gel phase and 
vice versa with sizes smaller than 0.01 /im 2 have been 
reported in a recent AFM study [l^, and have been re- 



4 



lated to lateral strain resulting from density differences 
in both phases. 

In summary, the high resolution neutron measure- 
ments show that the dip in the dispersion relation is 
actually deeper in the gel than in the fluid phase, ac- 
cording to the assignment of the peaks in e.g. Fig[21 (b) 
as based on their T-dependence. In this case, the present 
data also has important implications for the phase tran- 
sition from the gel to fluid phase and respective phase 
coexistence in lipid bilaycrs. Finally, we argue that the 
minimum £!o must be smaller in the gel phase, since in 
the framework of the GTEE model it can be shown that 
the frequency of the excitations is approximately equal 
to il — Q r vo/ y/2S(Q r ) (valid in the region of the mini- 
mum, see where vg = \JkBTjm s is the thermal 
velocity of the scatterers |22j. If this relation holds, it is 
clear that the minimum is deeper for the gel phase, where 
the elastic scattering S(Q r ) is more strongly peaked. Let 
us illustrate the relevant length and time scales of the 
excitations. With the given values of fin and the cor- 
responding wave vector Qo, we can estimate the wave 
velocity of the corresponding propagating modes in the 
two phases: cl/b = Qo/Qo — 100 m/s and clo. ~ 220 
m/s. From the measured width of the excitation, we es- 
timate the lifetime r = 2tt/Aui s ~ 17 ps and 9.2 ps, as 
well as the dynamic correlation length of the excitations 
£ T = ct ~ 17 A and 20 A, for the gel and fluid phase, 
respectively. These length scales are almost identical for 
the two phases, and fall very close to the static correla- 
tion length £ r!( 3 ~ 20A in the gel phase. 

It can be speculated that the dispersion relation in 
both phases is important for the transport of molecules 
across the bilayer as well as parallel to the bilayer. In 
particular the self-diffusion of the lipid molecules may be 
coupled to the collective motions of the lipid tails. Free 
area needed for diffusion might be generated by density 
waves which are incommensurate with the mean inter- 
molecular spacing, 2n/Q . The frequency Q/(2tt) of the 
collective density waves may serve as trial frequencies 
for thermally activated chain conformations which are 
involved in the formation of free volume. The free vol- 
ume model of diffusion gives a diffusion constant D = 
D'exp(— a*/af — Ea/iksT)) , which depends on the crit- 
ical free area a* per molecule needed for a hopping pro- 
cess, the average free area per molecule a/ and the en- 
ergy barrier E a for diffusion |l7tll8l |. Possibly, trial jump 
rates invoked in the model (in the pre-exponential term) 
may be related to the frequency of collective motions. 
Furthermore, the energy of the collective chain motions 
measured here, must be related to the average amplitude 



of excitations at thermal equilibrium, these amplitudes 
may determine the probability for the generation of crit- 
ical free volume. Certainly, more work both experimental 
and theoretical is needed to investigate the relationship 
between self-diffusion and collective chain motions. To 
this end, studying elementary diffusion events by molec- 
ular dynamics simulation may provide important insight. 
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